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SYNOPSIS 

ELECTRONIC STRUCTURE OF DISORDERED ALLOYS . 
USING MUFFIN-TIN POTENTIALS 
By 

S. S. RAJPUT 
Department of Physics, 

Indian Institute of Technology, Kanpur. India 


Disordered alloys lack the translational invariance and 
therefore, the electronic structure of these systems can not be 
studied by using standard methods of the band theory. However, 
a great deal of progress in this regard has been made using the 
coherent potential approximation (CPA> to the realistic 
muffin-tin potentials. The CPA is a mean field approximation in 
which the disordered alloy is replaced by an effective medium 
which IS determined by a self-consistent condition. This theory 
when applied to the muffin-tin model of disordered alloys. is 
known as Korr i nga-Kohn-Rostoker coherent-potential- 
approx imat i on <KKR-CPA> . This theory reduces to the standard 
Korr i nga-Kohn-Rostoker (KKR> band theory in the perfect crystal 
limit. The KKR-CPA offers a parameter-free first-principles 
calculation of the electronic structure of substitutional 
disordered alloys within the local density approx imat-i on of the 
density functional theory. 



X 


In chapter 2, a brief discussion of the KKR-CPA theory 
using the multiple scattering forrralism is presented. The 
solution of the multiple scattering equations is greatly 
facilitated by the muffin-tin approximation of the constituent 
potentials. The charge self-consistent KKR-CPA is then derived 
using the Green’s function approach. This is then applied to 
study the electronic structure of Nb Mo alloys for various 

X 1-x 

concentrations <x = 0,25. 0.50 and 0.75). Results of the 
self-consistent potentials, charge densities and density of 
states are presented and studied as a function of Nb 
concentration in the alloy. It is seen that the calculation of 
the density of states is very sensitive to slight changes in 
the potentials. The density of states at Fermi energy as a 
function of Mo concentration has been compared with 
experimental data obtained by the analysis of specific heat 
data and a good overall agreement is found between the two. The 
results of the density of states for Nb ^Mo ^ have been 
compared with soft x-ray spectra <SXS> curves and a very good 
agreement has been found between the experimental and fully 
charge self-consistent KKR-CPA results. Our results of the 
density of states for Nb ^Mo ^ have also been compared with 
earlier non-charge self-consistent KKR-CPA calculations and are 
found to be quite different in certain energy range. This 
underlines the importance of full charge self-consistency. We 
have also found some deviations from Nb-based rigid band model 



xi 

of Nb-Mo alloy. These deviations decrease with increasing Nb 
concentrat i on . 

The KKR-CPA is a single-site approximation and does not 
include the correlated scattering from clusters of atoms. This 
may be important in systems having clustering tendencies and 
short-range order or in magnetic systems where magnetic moment 
of an atom is sensitive to the local environment. There have 
been many attempts to generalize the CPA to induce correlated 
scattering from clusters, but these have been mainly focused on 
simple tight-binding models. One aim of this work is to achieve 
this objective using the muffin-tin model. This has been done 
by using the augmented space formalism of Mookerjee, which 
guarantees the herglotz properties of the average Green’s 
function. In chapter the augmented space formalism is first 
applied to the t i ght-bi nd i ng Hamiltonian and tight-binding 
cluster CPA (TB-CCPA) equations are derived. The generalization 
of the tight-binding method then forms the basis for extension 
of this method to the KKR framework and KKR cluster CPA 
<KKR-CCPA> equations are derived. We have also solved the 
two-impurity problem in the tight-binding framework to explain 
the structure in the cluster CPA (CCPA) density of states. 

In chapter 4, we have applied the KKR-CCPA formulation to 
the s-phase shift semicircular model, which has an analogue in 
the tight-binding framework. A one-to-one correspondence has 
been established between the KKR-CCPA equations for this 


mode 1 
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and the TB-CCPA equations for the analogous semicircular model. 
This IS one of the main results of our work. The TB-CCPA 
formulation is then applied to a Bethe lattice model which in 
certain limit reduces to the tight-binding semicircular model. 
We found that the difference in the CPA and CCPA density of 
states IS appreciable only when number of nearest neighbours Z 
IS small and decreases as Z increases. Also in the CPA the 
minority band is smooth where as in the CCPA it gains 
structure. The structure in the CCPA density of states is seen 
at energies close to the impurity levels. This clearly 
indicates that the structure appears due to the correlated 
scattering from the clusters embedded in an effective medium. 
For a large value of Z, there is a little difference in the CPA 
and CCPA density of states for this model and s-phase shift 
semicircular model. This is expected because of their 
equivalence in this limit. 

In chapter 5. the main conclusions are summarized. The 
various shortcomings and limitations of the present work are 
discussed. Some suggestions for further work in this direction 


are also given. 



CHAPTER-1 


INTRODUCTION 

Nature has provided us with a rich variety of disordered 
alloys having diverse electrical, magnetic, transport, optical 
and superconducting properties, which make them very useful for 
technological applications. To understand these properties a 
knowledge of their electronic structure is essential. The 
disordered alloys lack the translational invariance which is a 
characteristic of ordered solids. As a result the Bloch's 
theorem, which greatly simplifies the electronic structure 
calculation of ordered solids, is inapplicable to these 
systems. In other words, the standard band theory methods 
developed for ordered solids can not be applied to disordered 
a 1 1 oys . 

Substitutional random binary alloys are the simplest kind 

of disordered systems and will be the main subject of this 

thesis. In such alloys there exists a lattice structure but 

each lattice point can be occupied by either of the constituent 

atoms. Examples of such systems are a-brass Cu-Zn, Cu-Ge , 

Cu-Pd. Cu-Ni , Cu-Au, Ag-Pd, Nb-Mo etc. Such an alloy will be 

denoted by A B , where x and y are concentrations of A and B 
X y 

atoms in the alloy, respectively. 


During the last two 

decades, a great 

progress 

has 

been 

made in understanding 

the e 1 ectron i c 

structure 

of 

the 
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disordered alloys, by the application of coherent potential 

approximation <CPA)^ The CPA is a mean-field approximation, 

which for the sake of configurational averaging, replaces a 

disordered alloy (A B > by an ordered solid of effective atoms. 

X y 


These effective atoms 

are 

de t ermi ned 

by 

using 

the 

self-consistent condition 

that 1 f 

an A<B> atom 

1 s 

embedded 

1 n 

this effective medium, the 

average 

scatter i ng 

x: 

i 

respect 

to 


the medium is zero i.e. 


X 


^CPA 

A 


+ 


y 


tCPA 

B 


0 


<1 . 1 ) 


where t is atomic scattering operator for an A or B atom 

A( B > 

embedded in the effective medium. A simpler approximation is 
average t-matnx approximation <ATA)^ in which the t-matnx 
corresponding to the effective atom is 


ATA 

t =xt4*yt <1Z> 

A ^ B 

where t is atomic scattering operator for an isolated A<B) 

AxB/ 

atom. Another simple approximation is virtual crystal 

4 

approximation (VCA> in which the potential corresponding to 
the effective atoms is assumed to be 


VCA 

V 


XV + 

A 


y 


<1 .5) 


where v is the potential of A<B> atom. This approximation 

A<B> 
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IS good only when the difference between the two constituent 
potentials is small. Note that ail these approximations (CPA, 
ATA, VCA) are single~site approx i ma t i ons ; i.e. they neglect 
correlated scattering from clusters of atoms. The CPA has been 
found to be the best single-^site approximation for calculation 
of electronic structure of disordered alloys. 

A frequently quoted model in alloy theory is the rigid 
band model^. In this model the potentials of the constituent 
atoms are assumed to be the same. However, Fermi energy 

V 

adjusted to give the required number of valence electrons <N ) 
per atom in the alloy as 

E 

n'' = J p(E) dE (1.4) 

-00 

where p(E) is the density of states (DOS) per atom of the host 

system. The rigid band model is a crude approximation and does 

not explain many experimental results. But because of its 

simplicity, it is frequently used as a first step to get an 

idea of the electronic structure of the alloy. 

A solid contains a large number of interacting electrons 
interacting with nuclear potentials. Therefore, the 
calculation of energy levels of such a system is a many-body 
problem, which can not be solved exactly. However, this problem 
can be solved to a great degree of accuracy by using the one 
electron approximation. In this approximation, the many-body 
problem is reduced to a problem of one electron moving in some 
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effective potential which is determined se 1 f -cons i stent 1 y . 

Muffin-tin approximation is further used to simplify the 

calculation of this potential. In this approximation, the 

potential is assumed to be spherically symmetric within a 

sphere of radius r and constant in the space between this 

m 

sphere and Wigner-Seitz cell as 


v<r) 


f'^A < B > 

Iv" 

A < B > 


(r> for r < r 


m 


for r > r 


m 


The IS a constant interstitial potential as^ 


<1.5) 


I 

V 

A<B> 


VS 


3 J "" v(r> r*" dr / <r 


3 

vs 


m 



< 1 . 6 ) 


where v<r) is a spherically symmetric potential which, it is 
hoped, will be slowly varying in the region between the muffin 
tin sphere and the Wigner-Seitz sphere of radius r . For a 
pure A<B) solid, v* defines muffin-tin zero, i.e. this 
constant is subtracted from the potential <1.5> making 
interstitial potential zero. For the alloy, this is calculated 
as 

V* = X + y V* . <1.7) 

ALLOV A B 


Note that muffin-tin potentials are non-overlapping. The 
muffin-tin approximation is a reasonably good approximation for 
the potential as can be seen by referring to the book by 
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Moruzzi et al . Using this approximation, they have calculated 
various properties such as cohesive energy, bulk modulus. 
d";nsity of states and Fermi energy. They got very good 
agreements with experimental results and were able to explain 
general trends. Recent experience with alloys shows that the 

muffin-tin approximation is also a good approximation for 

. ,, ,, l-J.8-14 

meta 1 1 1 c a 1 1 oy s 

The muffin-tin potentials for the constituent atoms are 

constructed by using the local density approximation <LDA> of 

the density functional theory^^. The potential can be written 

5 

as 


V < r ) 

A< B > 


A{ B > . 

V < r > 

COUL 


V ( r ) + 

HUC 


A<B > . 

V < r ) 

xc 


< 1 . 8 ) 


where Coulomb, nuclear and 
are 


V 


A< B > 

COUL 


<r) 



exchange - corre 1 at i on 


p < r ’ ) 

A<B> 



contr 1 but i ons 


(1.9) 


v*‘®’ <?> = - e^ z /r . <1.10) 

HUC A<B> 


A< B 

V 

XC 


(r) = 


6E Ip (r>3 

xc *^A(B > 


6p <r) 

'^A<B> 


( 1 . 11 ) 


where p is the electron charge density in A<B> cell, z is 

Al { B 

the atomic number and n denotes the integral over the 
Wigner-Seitz unit cell. Here, is the exchange-correlation 

energy functional which in the local density approximation is 
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1 5 

given by 

^ rp<r> e <p(r>> dr , (1.12) 

where is the contribution of exchange and correlation to 

the total energy (per electron) in a homogeneous but 
interacting electron gas of density p. For paramagnetic case, 
von Barth-Hedin form of e 

xc 

(r > = ^ ^ 0.045 F<r /21) . <1.15> 

XC ® r m 

m 

where r^ is given by 

-r=— (r> = (1/p) 

y ® 

and function 

F(x> = (1 + x^> ln(l + -i-> - x^ + 

X 

Gunnarsson and Lundqvist form of c is^^ 

xc 

<r ) = — 0.035? FCr /11.4) . (1.16) 

XC or m 


1 8 

Vosko et al . form of e for paramagnetic case is 

X!C 



<1.14) 

- -i-. <1 . 15) 


where 



7 


2 7 \ / 7 1/2 

X(x) = X + b X + c. Q = (4c - b > and x = <r > ^ . (3.17) 

s 

Here, A =0.0621814. x = - 0.10498. b =5.72744 and c =12.9552. 

o 

It has been seen that use of different 

exchange-correlation potentials does not give large differences 

in the calculation of various properties'^. For example 

differences in cohesive energies in Li , Na , K and Rb were found 

18 7 

to be less than 8 mRy . Because Moruzzi et al . have found 
good agreement between theory and experiment for several metals 
including Nb and Mo. using von Barth-Hedin form of the 
exchange-correlation potential, we have used this form in our 
calculations as will be discussed in chapter 2. 

In the last decade. the CPA theory has been very 
successful as a single-site approximation to calculate the 
electronic structure of random substitutional disordered 
alloys. This theory when applied to muffin-tin model of 
disordered alloys. is known as Korr i nga-Kohn-Rostoker-CPA 
<KKR-CPA> theory^. In the perfect crystal limit, the KKR-CPA 
theory reduces to the standard KKR band theory^. In earlier 
stage of application, potentials v^ and v^ were not determined 
se 1 f -cons 1 stent 1 y . In the most sophisticated application of the 
KKR-CPA. V and v are determined se 1 f-cons i stent 1 y within the 

A B 

framework of the local density approximation described in 
chapter 2. This fully charge self-consistent KKR-CPA theory is 
a first-principles parameter-free theory of the electronic 
structure of random alloys. 
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The KKR-CPA is a single-site approximation and it does not 

include correlated scattering from clusters of atoms. This 

could play an important role in systems having short-range 

19-21 

order and clustering tendencies such as Cu-Pt and Cu-Pd. 

There have been many attempts to generalize the CPA to include 

22-24 

the correlated scattering from clusters of atoms but these 

have been mainly focused on simple tight-binding models. One 

aim of this work is to achieve this objective using muffin-tin 

potentials. This have been done by using the augmented space 

23 

formalism of Mookerjee which preserves the herglotz 

properties of Green's function. The augmented space formalism 
has been successfully applied to generalize the CPA theory in 
the tight-binding framework. Therefore, we first apply the 
augmented space formalism to the tight-binding Hamiltonian and 
tight-binding cluster CPA <TB-CCPA) equations are derived. For 
deriving the KKR cluster CPA (KKR-CCPA> equations we have shown 
that the Green's function (G .> in the tight-binding framework 
and the path operator <T } in the KKR framework are 

vj 

structurally similar. This similarity shows a one-to-one 
correspondence between the variables of tight-binding framework 
and the variables of KKR framework as shown in section 3.6. By 
using this correspondence we have derived the KKR-CCPA 
equat i ons . 

Implementation of the KKR-CCPA to realistic system is 
difficult due to its complexity and involvement of lengthy 
computation. So as a test case, we have applied the KKR-CCPA 
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using two-atom clusters to s-phase shift semicircular model. 

This model has the following advantages (i) The KKR-CPA theory 

2 5 

has been applied to this model in 1970 by Soven (ii> The 
model has an analogue in the tight-binding framework 
( t 1 ght-b 1 nd 1 ng semicircular model). A one-to-one correspondence 
has been established between the KKR-CCPA equations for s phase 
shift semicircular model and TB-CCPA equations for the 
analogous tight-binding semicircular model. This is one of the 
main results of our work. The tight-binding semicircular model 
is a limiting case of the Bethe lattice model, for the limit of 
large number of nearest neighbours. Therefore. the TB-CCPA 
formulation is applied to the Bethe lattice model and the 
structure of the CCPA density of states is explained on the 
basis of clusters of two atoms embedded in the effective 
med 1 urn. 

The outline of the thesis is as follows. In chapter 2, we 

2 6 

apply the charge self-consistent KKR-CPA to Nb-Mo alloys 
Section 2.1 presents the introduction to this chapter. In 
section 2.2, we derive the KKR-CPA equation in multiple 
scattering formalism. In section 2.5, we present the charge 
self-consistent KKR-CPA method. The Green’s function method is 
used to calculate charge densities of constituent atoms in the 
alloy. In section 2.4. we present computational details of our 
calculations. We have chosen von Barth-Hedin exchange 
correlation potential to calculate the new crystal potentials 
of constituent atoms in the alloy. Flow chart of the charge 
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se i f-cons 1 stent KKR-CPA is also presented in this section. The 

mixing schemes for the CPA and charge se 1 f- cons i stency loop are 

discussed. In section 2.5, we study the charge self-consistent 

potentials, charge densities and the density of states as a 

function of Nb concentration in the Nb Mo alloy for x=0.25, 

X y 

0.50 and 0.75. We find that maximum differences in the Nb<Mo) 
valence charge density and potential with respect to the pure 
Nb<Mo> valence charge density and potential, are of the order 
of 10^ and 5% respectively, near the muffin-tin radius. We note 
that these differences decrease with increasing concentration 
of Nb<Mo> in the alloy. In this section, we also compare the 
results of the charge self-consistent KKR-CPA density of states 
with the density of states using the Nb-based rigid band model 
and find differences in some energy range. For Nb ^Mo _ alloy, 
we compare the results of component density of states 
calculated by using the charge self-consistent KKR-CPA method 
with the results of soft x-ray spectra curves. We find that 
there is a good agreement between our results and experimental 
results. We also compare our results with the earlier 
non-charge self-consistent KKR-CPA calculations. In the final 
section 2.6, we summarize our main conclusions. 

In chapter 5. we present formulation of the CCPA. Section 
5.1 gives an introduction to this chapter. In section 5.2, we 
discuss briefly the augmented space formalism of Hooker jee. In 
section 5.5. the augmented space formalism has been applied to 
the tight-binding Hamiltonian and tight-binding CPA <TB-CPA) 
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equations are derived. In section 3.4, we derive the TB-CCPA 
equations by application of the augmented space formalism to 
the t 1 ght -b 1 nd 1 ng Hamiltonian. In section > . > . the two-impurity 
problem is solved in the tight-binding framework. In section 
3.6, the KKR-CCPA equations are derived with the help of the 
tight-binding CCPA <TB-CCPA) equations. Finally in section 3.7, 
the conclusions are presented. 

In chapter 4. we apply the KKR-CCPA and the TB-CCPA 
formulation to different models. Section 4.1 is an introduction 
to this chapter. In section 4.2. we discuss the s-phase shift 
semicircular model in the KKR framework. The path operators and 
formula for the KKR-CCPA density of states are derived in this 
section. In section 4.3, we briefly discuss the Bethe lattice 
model in the tight-binding framework. In section 4.4, the 
KKR-CCPA equations for the s-phase shift semicircular model and 
the TB-CCPA equations for the semicircular model are compared 
and a one-to-one correspondence between them is established 
<Table-4.1). In section 4.5, we apply the KKR-CPA and the 
TB-CCPA theories to the s-phase shift model and Bethe lattice 
model respectively. We also calculate the local density of 
states by embedding the two-atom cluster in the CPA medium. We 
give our conclusions in section 4.6. 

Finally, in chapter 5 we give the summary of achievements 
of the present work. Some suggestion for further improvements 
are given and some remaining problems are mentioned. 
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CHARGE SELF-CONSISTENT KKR-CPA AND ITS APPLICATION TO Nb-Mo ALLO' 

2.1 INTRODUCTION 

The Korr 1 nga-Kohn-Rostoker cluster coherent potential 

approximation <KKR-CPA) theory has been very successful in 

understanding the electronic structure of disordered alloys. As 

was mentioned in Chapter 1, it is basically an application of 

1-3 

the coherent potential approximation (CPA) to muffin-tin 

model of the disordered alloys. In the CPA. one replaces a 
disordered alloy by an ordered array of effective atoms which 
are determined by a self-consistent condition. The muffin-tin 
approximation of the potential plays an important role in 
simplifying the theory. Because the constituent muffin-tin 
potentials are spherically symmetric within a certain radius 
r . various quantities such as t-matrices can be expressed in 

m 

simple forms in the angular momentum space. The non-overlapping 

nature of the constituent potentials further simplifies the 

multiple scattering equations to tractable forms. 

The KKR-CPA reduces to standard KKR band theory in the 

perfect crystal limit. Within the local density approximation 

1 5 ■" 1 8 

of the density-functional theory (OFT), it can be made 

fully charge self-consistent. Thus charge self-consistent 
KKR-CPA provides a first-principles parameter- free theory of 
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electronic structure of disordered alloys. The theory rests on 
the same theoretical footing as the band theory of pure metals. 

The aim of this chapter is to study the electronic 
structure of disordered Nb-Mo alloys by the application of the 

p A 

charge se i f -cons i stent KKR-CPA . We have chosen Nb-Mo system 
because it is one of the important bcc refractory alloy 
systems. The system has relatively high superconducting 
transition temperature, high melting point and shows 

interesting behaviour in several properties such as electronic 

27 2 8 

specific heat and thermal expansion coefficient as a 

function of concentration. Also there are many experimental 
29-54 

results available on this system. To our knowledge, no 

fully charge self-consistent KKR-CPA calculation has been 
reported on this system till now. although some non-charge 
self-consistent KKR-CPA calculations have been reported^^ * . A 
very interesting behaviour of Nb ^Mo ^ was noted in these 
calculations. The density of states below Fermi energy showed 
approximately the rigid band behaviour, while the density of 
states above Fermi energy showed large deviation from the 
results of the rigid band model. This is surprising because 
this system is regarded as a classic virtual crystal system. 

The outline of this chapter is as follows. In section 2.2, 
the formulation of the KKR-CPA will be briefly discussed. In 
section 2.5, we present a brief formulation of the charge 
self-consistent KKR-CPA. In section 2.4, we discuss the 
application of this theory to Nb-Mo alloys. The computational 



14 


details are discussed in this section. In section 2.5, we have 
studied the charge self-consistent potentials, charge densities 
and the density of states as a function of Nb concentration in 
the alloy. Nb-Mo alloys with 25, 50 and 75^ of Nb are 
considered. We compare our results with the results obtained by 
using Nb-based rigid band model and earlier non-charge 
self-consistent KKR-CPA calculations. We also compare the 
density of states with the available experimental data and find 
good agreement between them. The final section 2.6, gives our 
main conclusions. 

2.2 KORRINGA-KOHN-ROSTOKER COHERENT-POTENTIAL-APPROXIMATION 

Consider a disordered substitutional alloy A B of A and B 

X y 

atoms with concentrations x and y respectively. The one 
electron potential in the alloy can be written as^ ^ 

N 

V(r) = E V (r > <2 . 1 > 

, ot % 

t = 1 

where the vectors r are defined by r. = r - ^ . 1^. are the 

t V V V 

location of the lattice points and N is the total number of 
sites in the solid. v^(r> <a=A. B> is the potential due to A or 
B atom and is assumed to be of the muffin-tin form and does not 
overlap with each other: 
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V (? > 

a \ 


V < r ) , 

a V 

constant , 


r < 

r > 



( 2 . 2 ) 


where r is the muffin~tin radius, 
m 

The Hamiltonian for a certain configuration of the alloy 
in atomic mass unit is 


H = - + V<r) . 


<2 . ?> 


The Green’s function in operator notation for a system of 
scatterers is defined as 


G = < El - H ) 


<2 .4) 


where E and I are energy and identity matrix. The equation 
<2.4> can be expanded using the Dyson’s equation as 


G = G+ GTG, <2. 5) 

o o ' o 


where G is the free electron Green's function and 
o 


T = V + V G T (2.6> 

’ o 


IS the total scattering operator. T rnay be written as 



16 


T = E T 


(.Z.7> 


where the operators are called path operators and satisfy 
the following multiple scattering equation 


T =t6 +tG ET 

V J I t J ° ^ 3 


(2 . 8 > 


Here , 


t = V < 1 + G t ) 

t L O L 


(2.9) 


IS the t-matrix that describes the scattering from an isolated 
potential on the i-th site. Equation (2.9> in angular momentum 
space has on-the-energy-she 1 1 matrix elements as 


t'" (k) = - < « ) ^ exp(i6 ) sin(<5 > , (2.10) 

A< B > l i 


1/2 

where 1 is the angular momentum index, x = < E > and <5^^ are 
phase shifts. 

The free electron Green's function G in real space can be 


G <r.r’> = 
o 


-IX V }(x r.) h <x r.) Y (r ) Y <r' ) 6 

"I >LnLmT 


+ E Y (r ) )(xr) B (ae r’ > Y (r’ >, <2. 11) 

" L n I n nm I' m L" m 


= r ,r’=r’ 4*1^ and r.( r.) is the smaller 

ri n mm s. .> 


where r 
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(greater) of the variable r and r’ . Also L = (l.m) is a 

r» m 

composite index, Y^(x) is the real spherical harmonics and 
and h^ are the spherical Bessel and Hankel functions 
respectively. The matrix B is defined as 


B^^'(«> 

nm 


An 1 H 


E 




I 1 


L 

1 


c 

i 


Y (^ -^ > 

L m n 
1 


where 


C 

LL. ' L 

i 


X ht <x|^ I > I ( 1 - <S > . 

I ^ n m J nm 

^ /s 

r dn Y (x)Y <x)Y (x> 

J XL L ' L 


( 2 . 12 ) 


(2.15) 


and h!^ denotes the outgoing Hankel function. 

Starting from equation (2.8) and after some algebraic 
manipulations, on-the-energy-she 1 1 matrix elements of the path 
operator can finally be expressed as 


T^^'(«) 

nm 

= 

L i t"^<«) - B («) 

3-1 

. 

nm 

(2.14) 

The Fourier 

transform of equation 

(2.12) 


C B,-^(x> 3 

K Li Li' 

= 

(1/N) r expC-ii<-(^ - ^ )3 B^^'<«) 

n m nm 

nm 

(2.15) 

is related 

to 

the well-known 

KKR 

structure 

f unct 1 ons^^ . 

equations <2 

.10) 

and (2.15) are 

very 

1 mportant 

in the KKR 


theory; t<w> depends only on the potential while Bj*’<w) depends 
only on the lattice structure. These two matrices play a 
central role in the KKR-CPA theory. 
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Now the Green's function G can be written as^ 
G<E.?.?’> = Z C(E.? ) Z™(E.r > 

3 LL. ' ^ ^ ^ 


-6 EZ^CE.r ) j’^<E.? ) . 

nm " L n L n 


( 2 . 16 ) 


where r and r' are within n th and m th muffin-tin spheres and 
the wave functions z"‘”’\E,r > and j"'""^<E.r > are 

L n<w> JL r»<Tn> 

respectively the regular and irregular solutions of the 

di fferential equa t i on 


r + V (? > -E 1 ) = 0 . 

L n<m> n<Tn> J L n<m> 

For spherically symmetric potentials one can write 


<2 . 17a) 



ZA<B > 

V 

<E. r> 

= Y <r> z^'®^ <E.r> 






Li t 


<2 . 17b) 


jA<B> 

<E. r> 

^ ^ . ACB) , ^ 

=Y<r)/ <E,r> . 



. A(B> 

where 

and 

A<B) 

are radial wave 

f unct 1 ons 

and are 

normal i zed 

such 

that for r ^ r they loin 
m 

smoothly to 


A( B > 



ft ^ 



<«.r> 

= <« 

r) i -IX 

l^A<B>J 

(x r) 

(2.17c) 


A( B > 

<^1 

(x , r ) 

= i,(« 

r > , 




-ri_ 1 j: j. ■ AiB> , 58 

The regular wave function z can be written as 

Xif 


A{B> ,A(B>.^. A{B>, 

<r. E) 


<2.17d) 


A. ( a ) 

where (p^ (E> is an energy dependent renormalization factor 


independent of r such that 


X n 58 

r f or r ► 0 . 
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Equations (2.16) and (2.14) are exact within the 

muffin-tin approximation and can be used to obtain the Green's 

59 

function for a cluster of atoms . However, here we are 
interested in an infinite system and would like to obtain 
ensemble average of G(E,r,r’). The ensemble averages of 
equation (2.16) are related to the ensemble averages of the 
path operators which will be determined by invoking the CPA. 
The CPA condition that the average scattering from each site 
must be zero, can be expressed in terms of Green's function 
operators as 

x<G> fv<G> = G, (2.18) 

A B 

where < G > denotes the restricted site average of G when 

A<B> 

0~th site IS occupied by an A<B) atom. Equivalently this 
condition can be written in terms of path operators as 


x<T > fy<T 

r»n n = A r»ri n=:B nr> 


(2.19) 


where < T > and < T > are the restricted site 

nn n=A nn n=B 

averages of T where n th site is occupied by an A<B) atom and 

nn 

T^ IS the nn path operator for the CPA medium. After solving 

nn 

equation (2.19) . we get the KKR-CPA condition as 


xCD] +y[D3 =I 


( 2 . 20 ) 


where 


c -1 -1 -1 

CD ] = C I + ( t ^ - t ) 3. . 

'■ A<B> LL' OO A<B> C LL 


( 2.21 ) 


and c IS used to label quantities for the CPA medium. The path 
operator matrix is given by 
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^ ^ ift - B-* (>e> I'^l dk . <2.22) 

c 

where O is the unit cell volume. We note that T involves a 

oo 

complicated Brillouin zone integration. Though B^ <»t) is an 
off-diagonal matrix, for cubic symmetry the integral reduces to 
a diagonal matrix^ for 1< 2 . By simple manipulations, one can 
reduce the KKR-CPA equation (2.20) with the help of equation 
(2.21) to a computationally simpler form as 



LL' 



+ 



+ (t 


-1 



) T 


oo 



.(2.25) 


This equation can be solved only by a numerical iterative 

-1 -1 

method. The input for iterating this equation are t , t via 
equation (2.10) and B+(j«) via equation (2.15). 


2.3 CHARGE SELF-CONSISTENT KKR-CPA 


For full charge self-consistency, the charge densities in 

2 

A and B cells are needed. These can be obtained . once the 
restricted site averages of the Green's function 
<G(E,r,r')> are known, as 

A<B> 

E 

p (r) = -<l/n) f Im < G<E.r.r’)> dE . (2.24) 

^ACB> J A<B> 

-00 

where p denotes the charge density associated with an A(B) 

^A(B> 

cell. Equation (2,24> can be further rewritten as 
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(r) = -<l/n) <?)Z^‘“Nr)D T^] dE . <2.25) 

A<B> J I A<B> OOjLL. 

-00 L J 


For the charge self-consistent KKR-CPA, the KKR-CPA 
equation (2.2>) is first solved for given potentials of A and B 
atoms. The new charge densities and are calculated from 
equation <2. 25). New potentials are then calculated by using 
the local density approximation of the dens i ty- f unct i ona 1 
theory as has been described in Chapter 1. The new potentials 
are used to solve the KKR-CPA equation and again a new set of 
potentials is calculated. We iterate this process until the 
potentials and charge densities get converged. 

To calculate the density of states and component density of 

Ai B Z 

states, we define the matrices F and F as 


pACB > 
LL ' 


<2.26) 


where Q denotes the integral over the unit ceil. Then the 
component density of states for the A and B atoms can be 
expressed as 

p <E) = -<l/n) Im r < G<E.?.r'> dr 

^A{B> A(B> 

or 

p (E> = -d/n) Im tr f T° 1„ <2. 27) 

where tr denotes the trace in L space. Since equation (2.27) is 

diagonal in L space, it allows L decomposed symmetry components 

<s, p, t and e ) of p <E> in the alloy. The cubic symmetry 
2g g a<b> 
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reduces the number of distinct elements four i.e. s, p, t and 

2g 

average density of states for the alloy is calculated 
as 

p(E> = -(l/T7)ImtrrCxF'^D + yF®D >1*^ 1 . <2.28) 

[_ A B OO Jl,L 

2.4 COMPUTATIONAL DETAILS 

We have applied the charge self-consistent KKR-CPA method 

to the Nb Mo alloy for different concentrations of Nb <x=0.25. 

X y 

0.50 and 0.75). These alloys have bcc structure and their 
lattice constants are presented^^ in Table 2,1. This systems 
forms a continuous solid solution in the whole concentration 
range . 


TABLE 2.1 


Lattice constants of Nb Mo alloys. 

X y 


S.No. 

Nb concentration (x) 

(■y=i-xj 

Lattice constants (a . u . > 

1 

0.00 

5.940 

2 

0.25 

6.015 

5 

0.50 

6.054 

4 

0.75 

6.155 

5 

1 .00 

6.250 
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For the charge se i f -cons i stent KKR-CPA. we start with 
atomic charge densities p^^ ) for Nb and Mo, and 

calculate the potentials > from equation (1.8) as 

A B 

explained in chapter 1. The t ^(»c) and t ^(«) are then computed 

A B 

for the input potentials (v^^, v^^). The KKR-CPA equation is 

A B 

solved iteratively within a tolerance of 0.0001, to calculate 

- 1 c; 

the t . Note that in each KKR-CPA iteration, T has to be 
c oo 

calculated using equation (2.22). which involves a complicated 

Bnllouin zone integration. In the initial stages of 

implementation of the KKR-CPA theory, this was a major problem, 

which was eventually solved by using the special direction 
41 

technique . Although this technique is quite simple, it does 

not give as accurate results as tetrahedron technique which has 

42 

been developed by Kaprzyk et al . for disordered alloys. 
Since, higher accuracy is essential for full charge 
self-consistency, we have used the tetrahedron method using 858 
k points in the irreducible part of the Bnllouin zone. We have 
observed that it is difficult to achieve convergence if we use 
less number of k points in the tetrahedron. The new charge 
densities ^ 3^^® obtained by equation <2.25>. The new 

potentials ^ alloy are then calculated using 

equation (1.8) of chapter 1. Again, these new potentials are 
taken as input to repeat the whole process. until the 
difference of old and new charge densities is within a certain 
tolerance (.£ ). This tolerance has been defined as 

CHABOB 
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€ 


A ( B > 
aUAMOK 


I 


m , out 
<P 

A ( B > 


t n . 

- p > 

( B > 


r ^ dr 


<2 .29) 


This tolerance is taken to be 0.0001 . The flow chart of the 

charge self-consistent KKR-CPA has been given in Chart 2.1. 

For the new charge densities, we have to calculate Fermi 
energy, as can be seen from equation (2.25). The Fermi energy 
is calculated using Lloyd formula for integrated density of 
states. The Lloyd formula^ is written as 

r/<E) = N*<E) + ^-IrnTx In detj t'^- t’^j+y In det| t^^ - |j 


- Im E In det i - Bj | , <2. 50) 

k 

where N^<E> is the integrated density of states and N*(E) is 

the free electron integrated density of states. Though this 

formula is easy to handle, it sometimes gives unphysical jumps 

which can occur in In <logarithm> of a complex quantity, used 

5 8 

in the formula. Kaprzyk and Bansil have derived a generalized 
Lloyd formula which does not give such unphysical jump and is 
given below: 

N^(E) = <Im/Ti) tr^ ^ E In C G"^(E,k> + - In <G^) 

ic 

+ X In C(v^*<r^>> ^ G^] + y In C<v'"<r^)> ^ 3 


<2 . ) 



25 


FLOW CHART OF CHARGE SELF-CONSISTENT KKR-CPA 



CHART 2.1 
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where 

G <2.52) 

A<B> C ^ ’ 

Here, the KKR-CPA Green’s function G and the associated 

c 

c 

1 og- de r 1 vat 1 ve (3 are the se i f-cons i stent solutions of the 
CPA equation (2.18): 


X G + y G = G 

A B C 


where 


<2.55) 


1 

N 


r C G~^ <E,k ) + (3^ - {f 

“ o 


-1 


(2.54) 


^A<B> ^E> 


d , A{B > 

- 3 — Z, 

dr l 


* r - r 

m 


< 2 . 55 ) 


and is the log-derivative for the spherical Bessel function. 

A AfBi 

Note that v' (r> and z <r) have been already defined in 

V V 

equations <2. 17). 

Although equation (2.31) also contains the logarithm 
functions, but unphysical jumps of these functions cancel with 
each other. The integrated density of states at the Fermi 
energy gives the total number of valence electrons in the 
alloy. Therefore, the Lloyd formula permits an evaluation of 
the Fermi energy without requiring an explicit computation of 
the density of states. We have used formula of Kaprzyk and 
Bansil. equation <2.31), to calculate the Fermi energy. 
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The calculation of charge densities from equation (2.Z>) 

using integration along real energy is very slow and time 

consuming. To speed up the calculations, the complex energy 
4 3,44 

method has been used to calculate the charge densities and 

potentials. The k-space integration is much faster in the 
complex plane because the integrand becomes smoother when the 
energy is complex. 

To obtain full charge self-consistency it takes a large 

number of iterations to achieve convergence. Therefore, to 

facilitate both the CPA and charge self-consistency loops, we 

7 45 - 1 

have used mixing schemes * . In the CPA loop new t is taken 

c 




ITCPA-4*! CPA 


[ c t ;' )-] 


[ C t;' 3°“*]: 


< 2 . 56 > 


where a is the mixing parameter for the CPA and itcpa is the 
number of iterations. In the charge se 1 f-cons i stency loop, the 
new potential is taken as 


r 1 = a'' f 

L A<B> JlTCHARaE+1 A<B> ^ 


A(B> IXTGHAROE 


t <1 


V ^ r c>ut I 

-a > 1 v I 

A<B> (B> J 


(B> IITCHAR0E 


(2 . 37) 


where 

A< B > 


IS the mixing parameter for A<B) potential and 
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ITCHAROE IS the number of iteration in charge density loop. We 
have taken between 0.90 to 0.98 for all concentrations of 
Nb in the alloy. The CPA mixing a is taken as 0.05. 

CPA 

The CPU time on HP-9000/850 system for running single 
charge self-consistency loop, is about 10 hours. The required 
number of iterations for the convergence of the charge 
densities in the charge self-consistency loop is about 25 for 
ail concentrations of Nb in the Nb-Mo alloys. The final 
calculation of density of states and component density of 
states IS done along the real energy axis and takes about 
another 15 hours of CPU time on this machine. 

2.5 RESULTS AND DISCUSSION 

Figure 2.1<a> shows the radial parts of the total and 
valence charge densities of pure Nb in good agreement with the 
calculations of Moruzzi et al.^. Note that there is slight 
difference in the value of the lattice constant used in the two 
calculations. We get excellent agreement with their calculation 
if we take the same value of the lattice constant. Figures 
2.1(b>, (c> and (d> show changes in the valence charge 
densities of Nb in the alloy with respect to pure Nb valence 
charge density for x=0.25, 0.50 and 0.75 respectively. As 
expected the maximur difference <of the order of 10% near 
muffin-tin radius) occurs for x=0.25, which decreases with the 


increase in x. 
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Figure 2.2(a) shows the self-consistent muffin-tin 
potential of pure Nb used in our calculation. Figures 2.2<b), 
(c) and (d) show changes in the self-consistent Nb potentials 
in the alloy with respect to the pure Nb potential for x=0.2?. 
0.5 and 0.75. This difference is maximum for x= 0.25 and is of 
the order of 5% near the muffin-tin radius and it decreases 
with the increase in concentration x. 

Figure 2.5(a) shows the radial parts of the total and 
valence charge densities of pure Mo in good agreement with the 
calculation of Moruzzi et al.^. Figures 2.5(b), <c) and <d) 
show how the Mo valence charge density changes in the alloy as 
a function of Nb concentration. The maximum change which is of 
the order of 10^ near the muffin-tin radius is observed for 
x=0.75, as IS expected. This change decreases with increase in 
the concentration of Mo in the alloy. 

Figure 2.4(a) shows the muffin-tin potential of pure Mo. 
Figures 2.4(b), (c) and <d) show how the Mo potential in the 
alloy changes as a function of Nb concentration. The maximum 
chang. 'which is of the order of 5^ near the muffin-tin radius 
occurs for x=0.75. This change decreases with the increase in 
the concentration of Mo (y) in the alloy. Thus, we find that 
the maximum dif-'erence in the valence charge densities of Nb 
and Mo in the alloy with respect to the pure Nb and Mo values, 
are of the order of 10* near the muffin-tin radius. The maximum 
difference in the potentials of Nb and Mo in the alloy with 
respect to the potentials of pure Nb and Mo, are of the order 
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(a) Nb25MOj5 


Energy { Rydberg ) 

Figure-2.? The average density of states by using the charge 

self-consistent KKR-CPA (full curves) and the density of states 

obtained by Nb-based rigio band model (broken curves) for the 

Nb Mo alloys for <a) x=0.2?. <b) x=0.50 ana <c) x=0.7>. The 

X y 

energy is measurea with respect to the Fermi energy. /w—i-X) 
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of 5% near the muffin-tin radius. Also, we find that, these 
differences for Nb (or Mo) decrease with the increase in the 
concentrations of Nb (or Mo) in the alloy. We also find that 
charge transfer from Mo to Nb , is of the order of 1 %. 

Figures 2.5(a), (b) and (c) show the comparison of the 
average density of states (DOS) calculated using the charge 
self-consistent KKR-CPA, with those of Nb-based rigid band 
model for x=0.25, 0.50 and 0.75 respectively. We find that 
below the Fermi energy, there is no qualitative difference 
between the two results for all concentrations. However. there 
IS a shift in the charge self-consistent KKR-CPA density of 
states towards the lower energy region with respect to the 
Nb-based rigid band density of states. This shift of density of 
states decreases with the increase in the concentration of Nb 
in the alloy. In the region above the Fermi energy, there is a 
marked difference between the charge self-consistent KKR-CPA 
density of states and the Nb-based rigid band density of states 
for all concentrations. This has also been observed by Giuliano 
et al.^^ when they compared their non-charge self-consistent 
KKR-CPA calculation with the results of the rigid band model. 
However, the deviation from Nb-based rigid band density of 
states is less in our calculations. 

Table 2.2 gives the Fermi energies (E^> for various alloy 
compositions and alloy density of states (DOS) at the Fermi 
energies calculated within the charge self-consistent KKR-CPA 
(present calculations) and non-charge self-consistent 


KKR-CPA 
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methods * . We see that the Fermi energies and the density of 
states at m the two calculations are not very different. 
However, we shall see below that the full density of states 
curves show substantial difference in the two calculations. 


TABLE 1,1 


Density of states (DOS) at Fermi energy (E > for Nb Mo alloys. 

F X y 


to 

z 

o 

X 

Present calculations 

Results of ref 40 

E 

F 

(Ry) 

DOS at E 

F 

St/<Ry . atom) 

E 

F 

(Ry) 

DOS at E 

F 

St/ ( Ry . atom) 

1 

0 . 00 

0.8114 

8.9 

0.826 

10.7 

2 

0 , 25 

0.7650 

7. l' 

0.764 

8.6 

5 

O 

O 

0.7485 

7.8 

1 

0.749 

9.6 

4 

0.75 

0.7052 

11.4 

0.709 

14.8 

5 

1 .00 

0.6880 

19.2 

0.682 

1 

19.5 

1 


Figure 2.6 shows comparison of the density of states (DOS) at 
Fermi energies. The two theoretical curves are quite similar. 
These results are compared with experimental results 


obta 1 ned 
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Figure-2.6 Density of states at Fermi 

of Mo in the Nb Mo alloys: present 

X y 

results of reference <40> (broken 
results of reference (27) (solid curve 


energy vs concentration 
results (full curve). 


curve > 
wi th 


and experimental 
circles). 
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by the analysis of the specific heat data. Our theoretical 
results are in good agreement with the experimental results and 
show the correct trend as a function of No concentration. 

Figures 2.7<a) and <b) show comparison of the charge 
self-consistent KKR-CPA density of states with two different 
non-charge self-consistent KKR-CPA results^^*^^ for x=0.50. The 
results of references 55 and 36 are different due to their 
different choice of input potentials for Nb and Mo. There are 
large differences between our results and the previous 
calculations. Though the potential differences of the Nb and Mo 
in the alloy with respect to the pure Nb and Mo potentials are 
small, yet there is substantial difference in the results of 
charge self-consistent and non-charge self-consistent 
calculations. This implies that the density of states is very 
sensitive to the changes in the potentials. These facts 
underline the importance of full charge self-consistency. 

In figures 2.8<a> and <b> . we show the comparison of our 
calculated component density of states for Nb and Mo with soft 
x-ray (emission and absorption) spectra curves^^ and non-charge 
self-consistent results of Donato et al.^^ for x=0.50. Above 
the Fermi energy, our results are very close to the results of 
Donato et ai . . But there are some differences below the Fermi 
energy. Our results are in better agreement with the 
experimental results. 

Figures 2.9<a, b. c, d. e). 2.10<a, b. c, d. e> and 

2.11<a. b, c, d, e) show the comparisons of s. p. t , e and 

2g g 
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Figure-2.8 <a) The component density of states of Nb and <b) 
the component density of states of Mo in the Nb alloy: 

present calculation (full curve). earlier non-charge 
self-consistent calculation of reference 06) (dashed curve) 
and experimental results of reference (36) (dashed-dot curve). 
The energy is measured with respect to the Fermi energy. 
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curves). The s. p. t , e and total components are shown i 

29 9 


figures (a), (b). (c), (d) and (e) respectively. The energy 1 

measured with respect to the Fermi energy. 
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Nb component density of states in the alloy with pure Nb 
results, for the concentrations x=0.2>. 0.50 and 0.75 
respectively. All curves of Nb component density of states m 
the alloy shift towards the low energy regions. This shifting 
of curves decreases with the increase in the concentration x. 
This IS because the change in Nb potential in figure 2.2, also 
decreases with increasing Nb concentration in the alloy. 


F 1 gures 

2.i2<a. 

b . c . d . e > , 2 . 

15(a. 

b . c , d , 

e ) 

and 

2 . 1 4<a , b. c . 

d. e) 

show the comparisons of 

B. P. 

e 

9 

and 

component density of 

states of Mo in 

the alloy with 

pure 

Mo 

results. for 

the 

concentrations x 

= 0.25. 

0.50 and 

0.75 

respect i ve ly . 

All curves of component 

density of states of 

Mo 

in the alloy 

shift towards the high energy 

reg 1 ons . 

But 

thi s 

shifting of 

curves 

decreases with 

the 

1 ncrease 

1 n 

the 

concentrat i on 

of Mo 

<y) in the alloy. 

Th 1 s 

also reflects 

the 

same trend as 

observed in figure 2.4. 






2.6 CONCLUSIONS 


We have calculated the self-consistent charge densities, 
self-consistent potentials, density of states and component 
density of states of the Nb^Mo^ alloy for x=0.25. 0.50 and 0.75 
by using the charge self-consistent KKR-CPA method. It is found 
that the maximum difference in the Nb (Mo> valence charge 
densities in the alloy with respect to pure Nb <Mo) is of the 
order of 10% near the muffin-tin radius. Also. the maximum 
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difference in the Nb (Mo> potentials in the alloy with respect 
to pure Nb (Mo) is of the order of 5^ near the muffin-tin 
radius. The differences decrease with the increase in the 
concentration of Nb(Mo> in the alloy. The results of charge 
self-consistent KKR-CPA, show deviations from the density of 
states (DOS) of Nb-based rigid band model. in some energy 
ranges. It is also found that there are significant differences 
between the charge self-consistent density of states and the 
density of states reported by non-charge self-consistent 
KKR-CPA methods. This implies that the density of states is 
very sensitive to the Nb and Mo potentials. Our calculated 
density of states for the alloy at is in good agreement with 
the experimental result obtained from the analysis of specific 
heat data. There is a good agreement between our results 
(component density of states) and soft x-ray spectra curves for 
Nb ^Mo ^ alloy. 



CHAPTER 3 


FORMULATION OF KORRINGA-KOHN-ROSTOKER CLUSTER-COHERENT- 

POTENT I AL- APPROX I MAT I ON* 

>.l INTROC CTION 

During the last decade, the coherent potential 

approximation (CPA) has emerged as the most successful single*" 

site approximation for calculating electronic structure of 

random substitutional binary alloys, both in the empirical 

tight-binding and first principles Korr i nga-Kohn-Ros toker (KKR) 

frameworks. Despite its success, the CPA, being a single site 

approximation, does not include correlated scattering from 

neighbouring sites. This could play an important role in 

systems having short range order and clustering 

19-21 

tendencies . To study the effect of correlated scattering, 

A 6 

Gonis et al . proposed the idea of an embedded cluster method. 
In this method, a cluster consisting of a central site and its 
shell of nearest neighbours was embedded in an effective medium 
determined within the Korr i nga-Kohn-Rostoker 

coherent-potent lal-approximation <KKR-CPA> . However this method 
IS not fully self-consistent. The idea of self consistent 

This work IS already published in 3. Phys.: Condensed Matter 


2. 265> (1990). 
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clu-ter CPA (CCPA> was proposed by some workers within the 
framework of some new approaches. like travelling cluster 
approximation ^ and Augmented space f orma i i sm^ ^ which 
preserve the herglotz properties of the Green’s f unct i on"^^ ’ . 
The augmented space formalism has been successfully applied to 
the calculation of electronic properties of binary alloys by 
several workers in the tight-binding f ramework^^ ' . 

One can use augmented space formalism to go beyond the 
single site approximation within the conventional KKR 
method . The aim of this chapter is to present this 
formulation, which will be called Korr i nga-Kohn-Rostoker 
cluster-coherent-potentiai-approximation (KKR-CCPA). Since 
there is a structural similarity between the tight-binding 
Green's function and KKR path operator, we shall first discuss 
augmented space formalism in the tight-binding framework and 
then discuss generalization to KKR framework. In section 5.2. 
we start with a brief discussion of the augmented space 
formalism in the tight-binding framework. This is first applied 
to a tight-binding Hamiltonian and tight-binding CPA <TB-CPA) 
equations are derived in section 5.3. We derive the 
tight-binding CCPA (TB-CCPA> equations by application of the 
augmented space formalism to the t i ght-bi nd i ng Hamiltonian in 
section 5.4. Here, for simplicity, we limit our cluster size to 
two-atom only. In the tight-binding framework we have refined 
the technique such that various quantities appearing in the 


CCPA equations are calculated 


Till now these 



51 


quantities were being calculated by summing over an infinite 

49, >0 

series 

One of the features of the CPA is that in the low 

concentration limit it gives the features of a single 

impurity . In the same spirit, we want to examine whether CCPA 

density of states in the low concentration limit reduces to the 

two ““ 1 mpur 1 1 y local density of states. The formulation of the 

two "" 1 mpur 1 1 y problem is given in section 5-5- In some earlier 

49 50 

works on the CCPA * as well as in the present work, we note 
rich structure in the CCPA density of states. We also wish to 
understand the structure in the CCPA density of states in the 
light of the two impurities local density of states. These 
structure have been explained by calculating the two- impur i ty 
local density of states in the tight-binding framework in 
chapter 4. 

The derivation of the TB-CCPA equations forms the basis 
for extension of this method to the KKR framework. By using the 
correspondence between variables in the tight-binding and KKR 
frameworks, we derive the KKR-CCPA equations with the help of 
the TB-CCPA equations in section ?.6. Finally, conclusions are 
given in section >.7. 

5.2 THE AUGMENTED SPACE FORMALISM 

The augmented space formalism is a technique of 
configuration averaging where configuration fluctuations are 
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systematically taken into account. The augmented space 
formalism was first introduced by Mookerjee^^ and later Gray 

24 

and Kaplan provided a detailed exposition. It has been 
successfully used for averaging one particle Green's function 
and functions related to two-particle Green's function. We will 
discuss this formulation in the tight-binding framework. 

We consider the following tight-binding Hamiltonian with 
no off-diagonal disorder: 


H = + , (5.1a> 

where E is the site-energy corresponding to site i, and V 

'■i 

IS hopping integral. .?>< = lv><v|) and T ( = |v><j|> are 
respectively projection and transfer operators in the Hilbert 
space X spanned by the site labelled basis Clv>>. The elements 
of Green's function for this system are 


G, = <t| <E I- H> ^ ij> . 

^ J 


(5.1b) 


For a random binary alloy A B . the random site energy E. 

X y t 

can be written in terms of a random occupation parameters n^ 
as , 


where 


E 


E n + E (1-n ) 

A L B 


ri if A 
t.0 if :.= B 


(5 . lc> 
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with a probability distribution 

P<n) = X 6(n-I)+ y 6(n ) 

wh 1 ch sat 1 sf 1 es 

P ( n ) > 0 

V V 

fp ( n > dn = 1 . 

We may express p (n > 

1. V 

P < n ) = - ( 1 /rr ) Im < | < n I - M ) "‘ | f° > 

^ V t V t V 

where M is an operator in the configuration space B of 
rank 2 spanned by I and with a representation 

= X + y 4-(x y> ( T + T ) . <5. 2) 

»• V V V V 

2 5 

The augmented space theorem then tells us that 
configuration average of a function is 

</ CCn }] > = < >. <5.3> 

t V 

where \4>> ~ fl |f°>, is a member of the basis = FI | f >, 

V *■ i- 

< s = 0.1 and p =1,2 ...2*^ > belonging to the configuration 

\ 

space J = n 0 of rank 2^. The function f- [fKll is an operator 

I, 

N 

in the augmented space v' = S<f ® j of rank Nx2 . 
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By the augmented space theorem, the eiements of 
configuration averaged Green’s function are 

< . (3.4a) 

The augmented space Hamiltonian H. is given by, 

H = E «> 4»' + 6EJ]5>®M+rrv T » ^ , (3.4b) 

where 6E = E^-E^ and S' is the identity operator in the 

configuration space. Note that equation (3.4a) is exact, but 
cannot be evaluated because of the large dimensionality of the 
augmented space. Therefore, we look for some approximation 
which will reduce the size of the augmented space. For this 
purpose, we partition the augmented space into a subspace I 
spanned by |^,^> and remaining subspace denoted by 1 1 . I f we 
take a one-atom cluster, we will get the CPA equation as 
derived in section 5.3. For the CCPA, we have taken a two-atom 
cluster, as presented in section 5,4. 

3.5 TIGHT-BINDING COHERENT-POTENTIAL-APPROXIMATION <TB-CPA> 

In this case, we take the one-atom cluster IS consisting of 

23 

the 0-th site. The space I is spanned by two basis vectors 
|0,f°> and |0.f*> . 


Now we partition H as 
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H = 


where 



H 

I 


H* 1 


H’ 




<3.5a> 




< 5 . 5b> 


where 


E = 


xE +yE,E=xE+yE and W = (x y >2 < E - E ) 

A B B A A B 


By the partition theorem. <G > may be written as the 

'•j 

representation of the resolvent of an operator H in the 
subspace I where 


H = H + H’(E ^ - H ) . 

I II 

With the help of equation <5.?). we get 

H = H + A < > . 

o o 


where 
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A = E E V„. G° V, 




Oj J k ko 


( 5 . 6 ) 


Since we need only the < 0 , f j j 
'' -1 

(E>5‘’ - H> .we again partition H as 


,0 > 

o 


e 1 ement 


where 


H = 


H 


^ 12 


H 


12 


H 


H = <E + A> 5^ 

1 of 

o 


<3 .7a) 


H = W T^o^i 
12 off 

O O 


<3 .7b> 


H = <E + A ) ^,1^.1 
2 f f 

o o 


<3 .7c> 


By the partition theorem 

<0.f^|<Ei^ - H)"‘|f^.0> = <01CEv5f - ^ H^^3"^10>. 

<3. 8> 

Let the coherent-potential Hamiltonian be 


H = ro'^+EE'^ 

• f f ^ O I a: ; J ’■ J 

V J 


(3.9) 


of 


then the Green’s function for the coherent-potential medium 
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<G> = <(E«?-H >~^>. 




T h © rna "trix ©ifimBnt. ^^0 ^ G ^ is given tDy 


= (E-o'-A)~^ (5.11) 

00 o 

and 

G = ( E I - H )~^ = ( E I -E - A . 

2 2 


Putting the value of A from equation we get 


G = < g“^ - E + <y ) . <3. 12) 

2 00 o 


By comparing equations (3.8) and <5.9>. and eliminating 
with the help of equation <3. 12) we get self consistent CPA 
equation in terms of the self energy: 


o'=E+W^ . <3.13) 

o — 

G~^ + O' - E 
oo o 

This IS identical to the CPA equations generated from 

more conventional methods. Note that this confirms the 

uniqueness of the CPA <in a single site case) as derived from 

various approaches such as mui t i pi e^scatter i ng^ embedding 

1 .54 


approaches etc. 



58 


5.4 TIGHT-BINDING CLUSTER- COHERENT- POTENTIAL APPROXIMATION (TB-CCPA) 


We shall take, as an example, the cluster ^ consisting of 

sites 0 and 1; i.e. we will include correlated scattering from 

two sites only. 

2 

The subspacG I is spanned by eight basis vectors |0,<^ >, 

11. >. |1.0 >. 11. 0> and ll.<^„. 

* OjL O Ox 

where 4> = 1 f° >. <P, = 1 f ° >, <P =lf‘ f° > and 4> = 1 f‘ f^>. 

O 1 ' O 1 ' O i ^oi o ± 

Now we partition H as 



C 5 . 14a> 


We now replace the Hamiltonians and H' by translat lonal ly 
symmetric effective Hamiltonians and H' . Their diagonal 
and off-diagonal elements are a and a respectively. Thus 

O 
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L “ tP-o., ^ v.,;-o.i vj J I 4>^ 4>^ <t>^^ 


L jP^O,l "■> J ^ ^1 ^oi 


<3 


E=xE+yE. E=xE_+yE and W = (x y) < E - E 


B A 


A B 


By the partition theorem, <G > may be written as 

vj 

representation of the resolvent of an operator H 
subspace I where 


H = H +H'(Es^-H )”‘h''*‘ 

I II 


With the help of equations <3. 14). we get 


H = H +1? < +JP ) + (? T r >l®r 5>. +5>. +^. +5>. 1 

I Poo O 1 ^oi 01 ^lO lO I ^ <P <p <p 

L J ^ ^O 1 Ol 

where 


- i.-. « „«f f<o . 1 > 

f = F F O' C a 

oo \b^ , Oj Jk 

J , ki«0 , 1 


ko 


r v r o ' * V 

^oi ■ OJ Jk 

J , k?=0 , 1 


ki 


and 


f = F F O' G , <y. 

^lO ^ " ij jk 

J , kp^o , 1 


ko 


.14b) 

> . 

the 
in the 


(3.15) 
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Since we need only the \<p,}> element of (E ^ - H) 

we again partition H as 


where 


H 

1 


H = 


H 


H 

•- 1 2 


H 


1 2 


H 


2 J 


E+? 


oo 


V 

oi oi 


V 

lO ^lO 


E+? 


oo 


H = 


12 



0 

0 


0 

0 


0 0 W ' 
0 0 0 ^ 


<>.16a) 


The matrix H is of rank 6 and can be obtained from equations 
2 

<3.15> . By the partition theorem. 


<v.0|<E S' - H) = <vlCE S - H - H (E - H h’' 3~‘b>. 

1 12 2 12 

<5. i6b) 


Here , 


H <£ S 
12 


H ) 
2 


-1 


H' 


12 


r w Q w 

H Q W 

i<$ 


W Q W 1 

<si 

H Q W 

11 


where 


Q = C(E - H >"^3 

V j 2 V J 
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We define a trans lat i ona 1 1 y symmetric effective 
Hamiltonian H* such that its resolvent is the average Green’s 
function <G > : 

i-j 

H*"” = o'j, E + E E j j . <3.17) 


Here and are diagonal and off-diagonai elements of . 
By comparing equations <5.16b> and (3.17). we get 


O' = E + W Q W 
o <s<s 


a = V + W Q W 

Ol Ol 


a = V + W Q W . 

lO iO 


(3.18) 


Here 


Q = 

<S<5 <5 


Q = r"‘ r w r"^ t r'‘ w r'* t r"^ . 

i lO 2 a Ol 4 5 lO <5 


Q = r"‘ r r"‘ w r‘‘ r R"^ w r”‘ t r"^ . 

<J 1 <S Oi 9 4 lO a 2 Ol 1 


(5.19) 


wi th 


Ol 


V + f 

Ol Ol 


•y = V + f 

lO lO ' lO 


and 
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R 


, = E - E - ? 


oo 


R = E - E - ? - r R ^ <Y 

^ OO 01 1 lO 


R -E-E-? -WR^W 

3 OO 2 


R = E - E - ? - r r“‘ sr 

^ OO lO 3 Oi 


R =e-E-? -wr*w 


OO 


R^ = E - E - ? 


OO 


r r“‘ t 

01 5 lO 


<5.20) 


We have calculated ? . f and f analytically by using 

^00^01 lO 

the partition theorem on . The is given as 


,eff 


= < E I - H*”> 


ff , -1 


By partitioning this in 2x2 cluster subspace, we get 


oo 

G-^^1 

Ol 


f 

fE-o 

o 

-a 1 

Ol 





f 

G^” 

OO 


h. 

-a 

^ to 

E-o ^ 

O'* 


. ^o 

^oo- 



Equation <5.21) then gives. 
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oo 




Oi 




iO 



(3.22) 


Note that and as given by (3.22) are exact. 

49 50 

The earlier workers ’ used to calculate these quantities 
approximately from (5.15) by recursion method. 

3.5 TWO- IMPURITY PROBLEM 


The study of the two impurity-problem is important while 
examining the structures in the CCPA density of states. Let us 
embed two impurities A and B at the sites labelled as 0 and 1 
in the CPA effective medium. Then the Hamiltonian is 


h" 


O' 

o 


E ^ + E E ^ ^ 

h !*o ^ ± t , 1 


^ £ P ^ £ P 

AO B i 


+ V(r ) . 

oi lO 


(5.25) 


The Hamiltonian for effective medium is given by 


H 


f 


= a 


f£ f 


£ y> + <7 


► f f 


E E 




(5.24) 
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only for nearest neighbours and for no 

off-diagonal disorder. The effective medium may 

be pure solid or the CPA medium and the corresponding values of 

diagonal self energies o are E or ively . 

o a<b> o ^ 

The value of is equal to V for both cases. The impurity 

and the medium Green’s functions are 



(3.25a) 
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e 


®f f 

♦ f£ 

O’ 

“O' 

o 

1 

# f f 

r- 

i 

E“0 

o 


0 0 0 


1 

. <3 


T 

here G is the transpose of G 

By simple partitioning <5.23) into a space spanned 
and 1 and the rest, we get 




= < V'^ - G ^ 


re*'" 

G*” 

OO 

OI 



G 

G 

^ lO 

OO ■ 


ef f 


W 




(3 


After solving equation <>.27), we get 


G 


oo 


ft G*^^C< a - E )G*^^+ a ] 

' B OO Oi O B lO 


ft ft -[(O' - E >G*^^+a 3C(o -E )G*^‘'+a ] 

A B O A Oi O B iO 


<3 


ft G*”[<o - E >G*^^+ a ] 

A OO lO O A OI 


ii 


ft ft -[(o -E )G*''^4a 3C(o -E )G*"+a ] 

' A ' B O B lO O A Oi 


. 26 ) 


by 0 


27) 


28a) 


(3.28b) 
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f 


where is the Green s function for the effective medium, and 


!^ff 


( Wji 1 . 4 

O' - V ) G 

1 oo 


ff 


and 


r? = 1 + ( O' 


ft = I + { a 


E > G*^^+ 

A OO 1 


f f , , eff 

) G + ( a 

B OO 1 


V ) G 


.ff 

lO 


V ) G 


,eff 


oi 


(5.28c> 


The local density of states (LDOS) p on the i-th site is then 
obtained by 


p = -d/n) Im < G*^*^ ) 

t V t 


{ 


0 or 1 


} 


Equation (5.28) reduces in the single impurity limit for 
E =o' as (O' =V . a=0, ft =1) 

BO i, ' B 


ff 


where 


G = D G 

oo A OO 


D = C 1 t (a 

A O 


E ) G*" r' 

A OO 


5 . 6 KORRINGA-KOHN-ROSTOKER CLUSTER-COHERENT-POTENT I AL-APPROX I MAT I ON 


In the KKR framework, an alloy is modeled by a random 

array of muffin-tin potentials of A or B type. The on-the-sheil 

1-5 

matrix elements of single muffin-tin t-matrices are 


t'^‘®’(«) exp<i6^> sin(6^> 


(5.29) 
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1/2 

where x = (|E|) and 6^^ are phase shifts. The path 
matrices are given by 


operator 


T ([ C - B<x) ] 


(5.50) 


where C = t and B(x) is the matrix of real space structure 

function B^^(x). Here, we have suppressed the angular momentum 

indices. We note that B(x) depends only on the lattice 

structure and does not contain any disorder while C<h> has a 

binary distribution. For a random binary alloy the random 
variable C can be written in terms of the random parameter n^ 
as 

C = C n + C(l-n), 

V A t B X, 


where is given by equation (5.ic>. Also we define an 
operator D as 

D = T~‘ = C - B 

= rc5* - rrB t = c +6cE^n - eeb t , 

V 

and 

6C = C “ C . 

A B 

A comparison of equations (5.50) and (5.1b) shows that 
is structurally similar to and hence the augmented space 
formulation developed in section 5.4 for the 

tight-binding-CCPA. can straightforwardly be applied to derive 
the KKR-CCPA equations^^ The relationship between 
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tight binding and KKR formulations becomes clear if we note the 
following correspondence. 


G < » T 

E 4 >. C 

V 4 > -B <5.51 > 


E < > 0 


The procedure for generating the CCPA by partitioning of 
the augmented space is identical to that described in section 
5.4. By using equation (5.51), we get the KKR-CCPA equations as 

C = C - W'Q' W 
<s<s 

S = B + W'Q’ W’ . 

Ol Ol 

3 = B + W’Q' W’ . (5.52) 

lO iO i<S 

Here . 

S ,. = B + b , 

>•1 >-j '■J 


Q’ 



Q- = R’"V’ R'"^w'R'"V’ r:"‘w’r:~V’ r; ‘ . 

1 ID 2 3 Ol 4 5 lO O 




<5. 55) 



And 
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"^1 - B + F • 

oi oi ^ oi 


B + r ’ 

»o ^ lO ’ 


C - ^ ’ 
oo 


<y, p, ^ey, 

OI i lO 


W’R’~‘W' 

2 


cy, p. 

lO 3 OI 


W’R' W’ 

4 


c - - r* R’ ‘r’ 

oo OI 5 lO 


<5.54> 


where 


C=xC +yC,C=xC+yC and W 

A B B A 


(X ( C - C ) 

A B 


Also, 


f ’ =: C 

^ OO 


f r f f 1 “i f f "j ”1 

Ol ^ oo J J 


'S + 

OI 


[ c - c{c}“ <oT • 


■s + 

lO 


\ oo J lo [ oo Oi \ OO j lO J 


< 5 . 55 > 
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Now the average electronic density of states for the alloy 
can be calculated by the formula^ 

6 p(E) = Pgjf<E) - p^(E> = Im Tr|^<T j. <5.56) 

where P <E> is the free electron density of states, p <E> is 
the density of states for the effective medium. N is the total 
number of unit cells in the solid and B-» is the Fourier 
transform of . The configuration averages in equation (5.56) 
can be calculated using the augmented space formalism as shown 
in section 5.1. 

5.7 CONCLUSION 

We have derived the TB-CCPA equations by application of 
the augmented space formalism to the tight-binding Hamiltonian. 
We were also able to refine the TB-CCPA method such that 
various quantities, which were being calculated approximately 
so far, are calculated exactly . We have formulated the 
two-impurity problem to calculate the local density of states 
in the tight binding framework. We have derived the KKR-CCPA 
equations with the help of the TB-CCPA equations using the 
correspondence between variables in the tight-binding and KKR 


f rameworks . 



CHAPTER 4 


application of 

KORR I NGA-KOHN-ROSTOKER CLUSTER-COHERENT-POTENT I AL- APPROX I MAT I ON* 
4.1 INTRODUCTION 

Because of its complexity. implementation of the 

Kornnga-Kohn-Rostoker cl uster-coherent-potent lal -approx i mat i on 

(KKR-CCPA) in realistic systems is difficult and involves 

lengthy computation. To our knowledge no successful 

implementation has been carried out on any realistic system to 

this date. In this chapter, as a test case, we have applied the 

KKR-CCPA to an s-phase-sh i f t semicircular model. The model has 

several advantages: <i) Because of the semicircular modeling, 

the involved i<-space integration required to obtain the 

site-diagonal path operator, is bypassed. This reduces the 

computational effort and allows us to concentrate more on the 

effect of correlated scattering, which is the principal aim of 

this work. <ii) The Korr i nga-Kohn-Rostoker coherent-potential- 

approx imat 1 on (KKR-CPA) based on this model has been tried out 

2 5 

before by Soven . This allows us to compare the earlier work, 
based on conventional methods, with our augmented space 

*This work is already published in 3. Phys.: Condensed Matter 


2, 265? <1990). 
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generalization. However fho u .. r 

nuwever, the ad hoc assumption of a 

semicircular density of states for the structure factor, allows 

us only to obtain the site-diagonal path operator. It does not 

tell us how to obtain the off-diagonal elements, which are also 

necessary in the c 1 uster-coherent-potent i a 1 -approx i mat i on 

<CCPA). In order to circumvent this difficulty, we have first 

shown an exact analogy in the mathematical structure in the 

CCPA equations of the KKR method for this model and the 

1 1 ght -b 1 nd 1 ng method for the semicircular model. The path 

operator in the KKR method has its counterpart in the Green's 

function in the tight-binding framework. We have also shown 

through the analogy, that the semicircular model is essentially 

a Bethe lattice approximation. This allows us to model the 

off-diagonal elements of the path operator. It also gives us a 

clear insight into the approximations involved in the model. 

In section 4.2, we apply the KKR-CCPA formulation to the 
s-phase-sh 1 f t semicircular model. The Bethe lattice model in 
the tight-binding framework, is examined in section 4.5. The 
KKR-CCPA equations for s-phase-shi f t semicircular model and the 
tight-binding CCPA <TB-CCPA) equations for the semicircular 
model are compared and one-to-one correspondence between them 
is established in section 4.4. In section 4.5, we apply the 
KKR-CCPA and the TB-CCPA methods to the s-phase-sh i ft 
semicircular model and Bethe lattice model respectively. We 
also apply the formulation of two-impurity problem to calculate 
the local density of states, by embedding the two-atom cluster 
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in the coherent-potential-approximation <CPA) medium. In this 

way. we explain the structure of the CCPA density of states. 
The results of our calculations are presented in this section. 
Finally we give our conclusions in section 4.6. 

4.2 s-PHASE SHIFT SEMICIRCULAR MODEL 


We shall illustrate our KKR-CCPA formalism by applying it 

to a simple model. This will also bring us in contact with the 

single-site muffin-tin KKR-CPA work of Soven^^. As in that 

work, we shall assume that only the s-phase shift dominates, so 

that B, t and T become scalars in the angular momentum space. 

2 5 

For this model, like Soven , we shall assume that the 
phase-shift has resonance behaviour and the cotangent of phase 
shift IS given by 




^a(b ) 
^a(b ) 


E 


(4.1) 


where E . . i s 

A (b ) 

hal f -wi dth . 

In order 
the structure 


the resonant energy and is 

to Simplify the problem further, we 

^ a. 25 
f unct 1 on 


resonance 


assume that 


B-* 

A 


Bf 

q 

H 


(4.2) 


has no explicit dependence on energy. Here B^(x) is the Fourier 

transform of B . From equation (3.29) of chapter 3 we get 
'•j 
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C(»() = 


cot6^ + 1 = -ft c + 1 


(4.5) 


where c is given by equation (4.1). 


The path operator 


T 


-jq- E exp(-iq,^^ ) 


1 

< c - > 

q 


(4.4) 


T 

oo 


1 

» N 


< c + B-*) 
q 


■^J<c + b>‘/°<b>db. 


where /°(b> is a distribution function given by 


/°(b> = -4- E <5( b - B-^ ) . 

N q 

q 


(4.5) 


We choose 


/*^(b) as 


25 


Then, 


/"<b) = 4- < 1 
= 0 . 


|b| < 1 ■ 

|b| > 1. . 


oo 


oo 


where 


oo 


ZCc+(c -1) 


1/2 


1 . 


(4,6) 


(4.7) 


S 1 mi 1 ar 1 y , 
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T = 

o 1 


T 

lO 


1 




J/’Ch) 


c + b > * db . 


(4.8a) 


where distribution function /’(b) is defined as. 


/■(b) 


N 


exp(-iq.^ ) Si b 

-*■ 


- B- ) 


We choose /’(b) such that 
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(4.8b) 


/ ’ (b> 


- ( 1 - b^ )‘^^ b. lb! < 1 ■ 

0 . |bl > 1 


(4.9) 


where X is a constant. The validity of this choice will be 
discussed in section 4.4. Hence 


T = T = - -i- r (— XI -b^)‘"'*(c + b)"^b db 
oi io « J , rr 


where . 


os. 


X . ,2 

-a— <9 ) 

8 oo 



(4.10a) 


From equation (4.8b> we get. 


B = f/' (b) b db 

oi J 



(4.10b) 


Taking the Fourier transform of (4.2) we get. 


B 


oi 


B 


40 


» B 


oi 
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From equation (5.?>) of chapter 3 we get, 


= Z = b + b 

oi lO oi oi 


where , 


b = b - n ft 

Ol lO ^ Ol 


< 4 . 11 ) 


The Fourier transform of (4.11) gives us 


f b*^ 
q q q 


and 


b-^ <ae> = J] b expCiq.l^ > = 




Oj 


Oj 


b exp( 1 q ) 

Ol Ol 


The diagonal elements of path operator within the CCPA are 
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given as 


^7 E < C + B-* )"^ 

OO « N q 


1 + 


ft exp( 1 q )' 

^ Ol Ol 

< c + > 

q 


-1 


Expanding by Taylor expansion near 0 and retaining only 

first order term we get. 


f 

OO 


/? j, 1 ^ o 

04 r+i ( c + b > db 


j + — 21 /'(b> 

OO X «* _ J 




OO 




2 c + 


2 c - 


( c 


] 


or 


the 
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T*^^- _ 1 •ff 

I ^ ^ - Q 

OO ^oo ’ 

where 

C= 3„o * ^ ^ <‘.12> 

The o f f ““ d 1 agona i elements of path operator are given 


oi 


« N 


J] exp< - iq.^ ) <c 4- 
o± 


B- ) 
q 




ft exp( iq.^ > 

oi OA 

< c + B-* ) 
q 


=■ T 

oi 



c + b > 


db 


2 ft_ 


- T 


oi 


Oi 


X 


1 + 


< c - 1 > 


1/2 


Thus 


T.®ff i eff 

T = T = - g 

04 40 X 04 


where 


off 

^04 


= a 


Oi 


+ Z 


o± 




( c 


1 > 


4/2 


] • 


(4.15) 


Note that the form of the function retains its 

herglotz properties even after this additional approximation. 

Finally the change in density of states (6D0S> per atom is 
given by the equation (5.56) of chapter 5. which is 
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6p(E ) 


I m 
n 



d C 

d~r 


-> 


^ -b^> 

^ q q 

q 


“1 


N 


d ( B-*- + b*^ ) 
q q 

dE 

< 4 . 14 ) 



where 


/ ^ ^ \ 

\ ^oo d E / 


oo 


= u + 


( c > 

- -C - C ) ' > • 

A B 


<4.25) 


U = 


dC 

“dE 


-•ff 

, C = - X c + i 


Here eff is used to denote the effective medium, and 


C = X C + y C 


From equation <4. 15) 


U 


eff 


[ 


r ^ + y ' 

<1-0 

<1-0 f 

1 

1 1 

G" 

^ 2 X 

<c -c ) 1 

G 

J 

«" A B -J 


A B >■ 

A 

B "* 


(4.16) 


where. c = x c + y c 

A «* 


and c is given by equation <4.1 >. Also c is the cotangent 

A<®> 

of phase shift for effective medium within the KKR-CPA or the 
KKR-CCPA. By equation <4.14) 


« p(E> = - -i- Im t F.rst part t Second part 3 
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where 


First part = 

oo 


(4.17) 


Second part = 


— T’ < P*f f Q , . ~1 d ( B-* +b-*) 

E < C - -b^ > q q 


<4.18) 


Second part can be written as 


Second part = 


N 


E / ( B-» + b- ) 
Q q 


<4.19) 


where function /< B-* ) is defined as 

q 


d<B-*) 

/<B-») = <«C-l+B-*') — TF— 

q q dE 


<4,20) 


Now the Taylor senes expansion of /< B^ ) near b^ = 0 

<up to first order terms only) equation <4. 19) gives 


Second part = 


2N E 


I 


q 


<B-»+i> 

— 3__ - ft <B-++i) exp<iq.^ ) 


< c+B-> ) 

q 


oi q 


ft exp< iq > 
* oi oi 

< c + B-*’ > 

<t 


Oi. 


<4.21 ) 


Equation <4.21) can be written as 


Sj * =2 * ^ * =4 * S 


Second part 
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where , 


5, = — ^ J /*^(b>(c + b) ^b db = - 


1 2E 


2E 


2 2 1/2 
2c - 1 + 2c (c^-l> 


]• 


‘ r'*'V^<b)(c + b) ^db = 

-1 E L 


2 2E 


^ / 2 , , 1/2 

c + < c -1 > 


1 r+I 


r 


2 E J 2 

-1 < c + b 


^ ^oiT , 2 , , <6 c^- 4 c> 

6c - 1 + — 


4E 


, 2 ,,1/2 

< c - 1 ) 


]■ 


ft /’<b) db iX. ^ r . 9 - 

S, = - ^ _ = _i£‘ r 2 c * ' ^ ^ - 


2E J , , X L, n2 2E 

-1 ( c + b ) 


( c 


] ] ■ 


ft. 


1 ^ r ^ 

- 1 ' oi lo _ 2 


'5 2E 


J /'(b><c+b> db = j^2c -l + 2c(c -1) j 


<4. 22) 


Substituting the values of equations (4.16), (4.17) and (4.22) 

in equation (4.14), we get the formula for the change in 

52 

density of states as 


6p(E> 



(4.25) 


Note that this formula is valid for the KKR-CCPA as well 

• ff 

as the KKR-CPA depending upon the value of c. Here g^^ is 
given by equation (4.12). 
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4.3 BETHE LATTICE MODEL 


Bethe lattices have no closed loops and are completely 
characterized by their number of nearest neighbours Z or 
connectivity K = Z- 1 . For the tight-binding Hamiltonian on the 
Bethe lattice the Green’s functions for a perfect system can be 
calculated exactly by renormalized perturbation expansion^^ and 
are given by 


G == 
oo 


IK 

(K-l > (E-E^) + <K+l >■*'(£_£ >^-4KV^ 

o 


2V 

G = G = G r 1 . <4.24) 

I <E-Eo>+V^e-E >^-4K'/ 


where E and V are the diagonal and of f-di agonal elements of 
o 

the nearest neighbour Hamiltonian. We take V = O.^/VK which 
gives unit half band-width. For an alloy take the site 

energies E as E = *^5* and E = — 5- where S is band separation 

O At B C 

parameter. For the CPA medium, equation (4.24) reduces to. 


G = 


2K 


(K-lXE-o 


2V 


G = G = G 

0 1 lO OO 


o 




<4.25) 
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where cf is the diagonal element of the effective Hamiltonian 
o 

within CPA. For K >> 1 , equations <4.25) reduce to the Green’s 
function of the semicircular model as^^ 


G = ZTiE-a ^2 T 

oo 1 o (E-c? ) - 1 
L o 


G = = 

oi io 2 Tn< 


( G f . 

oo 


<4. 26) 


The Green's functions within the CCPA can be obtained from 
equatio. (4.25) by replacing V by where <7^ is the 
off-diagonal element of the effective Hamiltonian. 


4.4 COMPARISON OF KKR-CCPA AND TB-CCPA EQUATIONS 


In this section, we have compared the KKR-CCPA equations 
for the s-phase-sh 1 f t semicircular model of unit resonance half 
width with the TB-CCPA equations for the semicircular model of 
unit half band-width. Therefore, we have calculated the values 
of Green's functions for tight-binding semicircular model with 
the help of distribution functions. The Green's function within 

the CPA IS given^ as 


G. . 

'-i 



exp( - i ^ j ^ 



where 


<4,27) 
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‘^oi " ^lo = J b)"^ db , 

where 

/'(b) = — ^ J] exp(-iq,^ ) 6(b - S-+ ) . 

N 01 a 

q 

With the help of distribution functions /°(b> and /'(b) as 

defined in equations (4.6) and (4.9). we get diagonal and 

52 

off-diagonal elements of Green’s function as 

G = - r/®(b) (c'+ b)~‘ db = - 2C c'+ (c’*- . (4.28) 

00 

G = G = r/’(b) (c'+ b)"‘ db ^ ( G )^ . (4.29) 

01 iO j 0 OO 

where, 

c ' = O' - E 
o 
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Thiese equations for G and G are analogous to the 

OO Oi 

equations (4.26) which are exact and are derived by taking the 

large K limit of the Bethe lattice model. Since T in equation 

vj 

(4.4) and G in equation (4,27) are structurally similar, we 

'-j 

expect T and T in equations (4.7) and (4.10a) respectively 

OO Ol 

to have the form of equations (4.28) and (4.29). This is 
exactly what we get by using the distribution functions / (b) 
and /’(b). Therefore, we are justified in using these 
distribution functions. We also get 


V = V = f/ ’ (b) b db 

Ol lO J 


X 

8 


(4.30) 


The Green's function for the effective medium is given by 


>f f 


'■ J 


N 


exp(-lq.i^. > 

^ J 

E - O' -S-» 
o q 


(4.31a) 


where 


S-* - S-» + s-» 
q q 


c.* = V exp( 1 q .F^ . ^ • 

q Ol 


<4. 31b) 


>■ j 


y = O “V 

i j V J V J 


By using the 


distribution functions /“(b) and /’tb). «s hava 

— «ff 5 2 

off-diagonal elements of G as 


calculated the diagonal and 




By using these variables, the KKR-CCPA equations are 
compared with TB-CCPA equations in Table 4.1. It is clear that 
there is a one-to-one correspondence between them. 

4.5 RESULTS AND DISCUSSION 

We have applied the KKR-CPA and KKR-CCPA theory to the 
s-phase-shift semi c i rcu lar mode 1 , The change in density of 
states <6D0S> with respect to the free electron-gas. has been 

for this model. We have also 


calculated using these theories 
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COMPARISON BETWEEN KKR-CCPA AND TB-CCPA EQUATIONS 


KKR-CCPA EQUATIONS 
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calculated density of states within the t i ght -b i nd i ng 

coherent-potential-approximation <TB-CPA) and t i ght -b i nd i ng 
c 1 ust er-coherent-potent 1 al -approx mat 1 on (TB-CCPA> for Bethe 

lattice model for different K, 6 and concentration x. 

For convenience, we shall first present the calculations 
using the tight-binding method for the Bethe lattice model with 
K=2. Figure 4.1 shows the CPA and CCPA density of states for 
(5-1 in the low concentration regime. We note that there is not 
much difference in the ma)ority band between the CPA and CCPA. 
in contrast to considerable difference in the impurity band. 
Two sub-bands appear in the CCPA as compared to one in the CPA. 
These two sub-bands arise because of correlated scattering 
from the clusters of two impurities embedded in the effective 
medium. In the low concentration limit <x<<l>. the impurity 
pair in the alloy is most likely to be of AA type. The energy 
levels for the AA impurity pair embedded in pure B are at 
energies 0.55 and 1.04, as shown by the solid arrows in figure 
1. These two levels are close to the shoulder in CCPA density 
of states at about E = 0.6 and the peak at E = 1.0 

respectively. The peak at about E = 0.8 in the CCPA density of 
states arises due to the single A impurity level which is at E 
= 0.85 and about which the CPA impurity band is centered. Thus, 
the structures in the CCPA minority band density of states are 
accounted for by the examination of impurity levels. 

Figure 4.2<a) shows the CPA and CCPA density of states for 
<5 = 1 for a concentrated alloy (x = 0.5). We note that although CPA 




ENERGY (half band width} 


Figure ''.1 The averaged density of states using the CPA (dashed 
curve) and the CCPA (solid curve) for the parameters K“2. <5=1 
and X = 0.0>. Dashed and solid arrows indicate the positions of 
the single and two impurities levels in a pure B medium. The 
energy is in units of half-band width. 
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Figure 4.2 (a) The averaged density of states using the CPA 
(dashed curve) and the CCPA (solid curve) for parameters K=2, 
<S = 1 and x = 0.5. (b) Local density of states for two impurities 
pairs embedded in the CPA medium for the same parameters. Solid 
curves correspond to BA and AB pairs while the dashed curves 
correspond to BB and AA pairs. 





tjansity of states (arb units) 


ENERGY (half band width) 


Figure A.'} The averaged density of states using the CPA (dashed 
curve)anci the CCPA (solid curve) for parameters K = 2, 6= 1.7 ana 
x = 0 . ? . 
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density of states is smooth . the CCPA density of states has 
two peaks. We explain this on the basis of clusters of two 
atoms embedded in the effective medium. Unlike the low 
concentration limit, the probabilities of AA. AB . BA and BB 
pairs are identical . The peaks in CCPA around E = -1.0 and 
1.0 appear approximately at energies close to the local density 
of states peaks corresponding to BA and AB pairs embedded in 
the CPA medium as shown figure 4.2(b). The local density of 
states peaks corresponding to AA and BB pairs appear near the 
edges of CCPA density of states. This indicates that the extra 
peaks arising in the CCPA density of states are due to the 
correlated scattering from the impurity clusters in the 
effective medium. 

Figure 4.3 shows the density of states for <5=1,7 and 
x=0.5. We note that critical value of 6 for band separation in 
the CCPA IS 1.7 while in CPA it is 1.5. This implies that in 
the CCPA. the bands are less likely to be separated. 

The results for the Bethe lattice model with K=10 are 
shown in figures 4.4, 4.3<a). 4.5(b) and 4.6. In figure 4.4 we 
show the CPA and CCPA density of states for 6 = 1 in the dilute 


limit ( X 

= 0.05) . There 

IS no difference in 

the 

CPA and 

CCPA 

dens i ty 

of states 

i n 

majority band. The 

CCPA 

minor i ty 

band 

dens i ty 

of states 

has 

two peaks at about 

E = 0.6 

and 0.85 

1 n 


contrast to a very smooth CPA minority band. The two peaks in 
the CCPA density of states are well accounted for by the energy 
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Figure 4-4 The averaged density of states using the CPA (dashed 
curve) and the CCPA (solid curve) for parameters K=10. 6=1 ana 
x-0.05. The dashed and solid arrows indicate the positions of 
the single and two impurities levels in a pure B medium. 
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Figure 4.? <a) The averaged density of states using the CPA 
(dashed curve) and the CCPA (solid curve) for parameters K=10. 
6=1 and x=0.?. <b) Local density of states for two impurities 
pairs embedded in the CPA medium for the same parameters. Solid 
curves correspond to BA and AB pairs while the dashed curves 
correspond to BB and AA pairs. 
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levels of AA impurity pairs which are at E=0.64 and 0.87 
respect i ve 1 y . 

Figure 4.5(a) shows the density of states curves for <5=1 
in the concentrated limit (x=0.5>. We note that both the CPA as 
well as the CCPA density of states are smooth. This arises due 
to the following reasons : <i> The local density of states for 
the impurity pairs, as shown in figure 4.5(b), are smooth as 
compared to K=2. (ii> The local density of states peaks for 
(BA.AB) and (AB,AA) pairs are at about the same energies. 

Figure 4.6 shows the density of states for <5=1.4 and x 
=0.5. This again shows that the critical value of 6 in the CCPA 
is greater than that in the CPA. It is 1.4 in the CCPA where as 
it IS 1.2 in the CPA . 

A contrast between K=2 and K=10 cases is worth noting. In 
the dilute limit, the minority band for K=2 (figure 4.1) 
breaks up into two sub-bands while for K=10 (figure 4.4) we 
have only one impurity band with some structures. The 
separation between the two impurity levels for K=2 case is much 
larger than that for K=10. This difference is due to the 
different values of hopping parameter V = 0,5/yK, in the two 
cases. For K=10 case it is much smaller compared to that for 
K=2 . Thus in the concentrated limit, there is an appreciable 
difference in the CPA and CCPA density of states for K=2 
compared to K=10 case. The difference between the CPA and CCPA 
density of states decreases as K increases and eventually 
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reduces to zero when K is very large (=sl00). Note that in this 

limit, our model reduces to the semi -c i rcular model^^ as shown 

in section 4.5. 

Now we present the KKR-CCPA results for the s-phase-sh i f t 
semicircular model. We have calculated the change in density of 
states <6D05) with respect to the free electron gas within the 
CPA and CCPA. We have found that for the small X. there is 
negligible difference in the 6D0S between the CPA and CCPA, as 
was mentioned above for the Bethe lattice model with large K. 
This IB because for large K (or small X), CCPA equations for 
both models have one-to-one correspondence. But for large X 
there is a visible difference between the 6D0S in the CPA and 
CCPA as can be seen in figures 4.7(a) and 4.7(b>. Figure 4.7(a> 
shows the change in the density of states (<5D0S) using KKR-CPA 
and KKR-CCPA for same resonance half-widths (T =r =1) while 

A B 

figure 4.7(b) is for different half-widths (F =1, F =2). Note 

A B 

that these results, in general, are also valid for the 

tight-binding semicircular model as is evident from Table 4.1. 

As in the tight-binding case, to understand the structure 

in the KKR-CCPA density of states, we have to solve the two 

impurities problem. This requires a knowledge of the impurity 
2 

wave functions , which are not available for the s-phase shift 
semicircular model. However, since we have already shown a 
one-to-one correspondence between the KKR-CCPA and the TB-CCPA 
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equations, we expect that the mam features of the TB model 
will be reflected in the KKR model, and hence, the structure in 
figure 4.7 can be related to the correlated scattering from 
clusters of atoms. 

4.6 CONCLUSIONS 

We have applied the KKR-CCPA formulation of chapter 5 to 
the s-phase-sh 1 f t semicircular model, which has an analogue in 
the tight-binding framework. A one-to-one correspondence has 
been established between the KKR-CCPA equations for this model 
and the TB-CCPA equations for the analogous semicircular 
model. The TB-CCPA formulation (described in chapter 5> was 
then applied to the Bethe lattice model. We found that the 
difference in the CPA and CCPA density of states is appreciable 
only when the number of nearest neighbours Z is small and 
decreases as Z increases. Also, in the CPA, the minority band 
is smooth whereas in the CCPA it gains structure. The structure 
in the CCPA density of states is seen at energies close to the 
impurity levels. This clearly indicates that the structure 
appears due to the correlated scattering from the clusters 
embedded in an effective medium. For large value of Z, there is 
little difference in the CPA and the CCPA density of states for 
this model and the s-phase-sh i ft semicircular model. This is 
expected because of their equivalence in this limit. 


CHAPTER 5 


SUMMARY AND CONCLUSIONS 

In this thesis we have discussed and applied two methods 
for electronic structure calculation of disordered random 

alloys using muffin-tin potentials. The first method is the 
charge self-consistent KKR-CPA which treats disorder scattering 
in the single-site approximation and thus ignores the 
correlated scattering from clusters of atoms. The second method 
is the KKR-CCPA which is cluster generalization of the KKR-CPA 
method. The KKR-CCPA includes the correlated scattering from 
clusters of atoms, which may be important in systems showing 
clustering or ordering tendencies such as Cu-Ni, Cu-Pt and 

Cu-Pd or in magnetic alloys such as Ni-Cu and Fe-Ni in which 
the magnetic moment of an atom depends on the local environment 
of the atom. 

The charge self-consistent KKR-CPA was applied to Nb-Mo 
system which is a classic bcc refractory alloy system. We have 
calculated self-consistent charge densities, self-consistent 
potentials, component density of states and density of states 

of Nb Mo alloys for x=0.25, 0.50 and 0.75. It has been found 

X y 

that the maximum difference in the Nb<Mo) valence charge 
densities in the alloy with respect to the pure Nb(Mo> is of 
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the order of 10^ near the muffin-tin radius. Also, the maximum 
difference in the Nb(Mo) potentials in the alloy with respect 
to pure Nb(Mo) is of the order of near the muffin-tin 
radius. These differences decrease with increasing 
concentration of Nb(Mo) in the alloy. A comparison of the 
results of the charge self-consistent KKR-CPA and Nb-based 
rigid band model shows that the charge self-consistent KKR-CPA 
density of states deviates from the rigid band model density of 
states in some energy range. We also found that there is an 
appreciable difference between the density of states calculated 
by using the charge self-consistent KKR-CPA method and 
non-charge self-consistent KKR-CPA methods^ ^ . This shows 
that the density of states is very sensitive to the constituent 
potentials. We have found that charge self-consistent KKR-CPA 
density of states at Fermi energy different 
concentrations of alloys is in good agreement with the 
experimental results obtained by the analysis of specific heat 
data. We have also compared the results of component density of 
states calculated by using the charge self-consistent KKR-CPA 
method with the results of soft X-ray (emission and absorption) 
spectra for Nb _Mo _ alloy. We have found that there is a good 
agreement between our results and experimental results. These 
results underline the importance of charge self-consistency in 


the KKR-CPA theory. 
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The KKR-CPA calculation can ho 1tY^r^v^n _i 

idiion can be improved by going beyond 

the single-site approximation and using the KKR-CCPA which has 

been formulated by us. In this formulation we have used the 

augmented space formalism which guarantees the herglotz 

properties of the averaged Green’s function. The KKR-CCPA 

equations were derived bv first nntir%r^ ^ c ^ 

uy Tirsx noting a formal similarity 

between the multiple scattering equations of KKR theory and 

Green’s function of t i ght -bi ndi ng Hamiltonian, Therefore. we 

have first developed the tight-binding c luster-coherent- 

potent i a 1 -approx i mat i on <TB-CCPA) formulation. We were able to 

refine the TB-CCPA method such that various quantities which 

were hitherto being calculated approximately are calculated 

exact 1 y . In the tight-binding framework, we have formulated the 

two - 1 mpu r ity problem to calculate the local density of states 

on an impurity of two-impurity cluster. The KKR-CCPA equations 

were then derived by using the correspondence between the 

variables of KKR framework and tight-binding framework. 

The KKR-CCPA theory was applied to the s-phase shift 
semicircular model, which has its analogue in the tight-binding 
semicircular model. We have compared the KKR-CCPA equations for 
the s-phase shift semicircular model and the TB-CCPA equations 
for the t i ght-b 1 nd i ng semicircular model and a one-to-one 
correspondence between them has been established. We have also 
applied the TB-CCPA method to the Bethe lattice model. We have 
found that the difference in the TB-CPA density of states and 
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the TB”CCPA density of states is appreciable only when the 
number of nearest neighbours (Z) is small and decreases as 2 
increases. The minority band in the TB-CPA density of states is 
smooth and it gams structure in the TB-CCPA. The structure in 
the TB-CCPA density of states is seen at energies close to the 
two-impurity levels. This clearly indicates that the structure 
appears due to the correlated scattering from the clusters of 
atoms in the effective medium. For large Z. there is little 
difference in the coherent-potent lal -approx imat ion density of 
states and the c 1 uster-coherent-pot ent i al -approx imat i on density 
of states for this model as well as the s-phase shift 
semicircular model. This is because for large Z the Bethe 
lattice model reduces in the semicircular model as shown in 
section 4 . ? . 

In this work, we have applied the KKR-CCPA theory for 
two-atom clusters to a model case. The fully se i f -cons i stent 
KKR-CCPA theory has been developed by Razee et al.^^ in Green's 
function formalism. This theory can be applied to realistic 
alloy systems by taking 9-atom (bcc) and 15-atom (fee) 
clusters. The short-range-order can also be treated within this 
theory and charge self-consistency can be obtained. Such 
calculations are in progress in our group. 

The KKR-CPA calculation can be improved in many other 
ways. One should improve upon the local density approximation 
which we had used in our calculation. At present many workers 
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are engaged in this direction and have suggested various 

approximation such as gradient correction, GW approximation, 

weighted density <WD> approximation, and se 1 f - i nteract i on 

corrected (SIC) approximation etc.^^’^^. We feel that such 

approximations should be first tried and thoroughly tested on 

simple systems and pure metals before these could be applied to 

disordered alloys. The relativistic version of the KKR-CPA 

theory has been formulated by Staunton et al. and can be used 

to improve our KKR — results. However, such corrections in 

case of Nb-Mo alloys are supposed to be quite small. The other 

direction in which the theory for disordered alloys can be 

improved is by going beyond the muffin-tin approximation. The 

KKR-band theory for pure metals has been generalized to 

59 60 

non-muff 1 n-t 1 n potentials by Brown and Ciftan , Gonis and 
Faulkner^^. However, their methods have convergence problems as 
discussed by Gonis and Faulkner. To our knowledge the KKR-CPA 
theory has not been formulated for non-muff i n-t i n potentials so 
far and is still an open problem. 
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